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Broken discs: warp propagation in accretion discs 
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ABSTRACT 

We simulate the viscous evolution of an accretion disc around a spinning black hole. In general 
any such disc is misaligned, and warped by the Lense-Thirring effect. Unlike previous studies 
we use effective viscosities constrained to be consistent with the internal fluid dynamics of 
the disc. We find that nonlinear fluid effects, which reduce the effective viscosities in warped 
regions, can promote the breaking of the disc into two distinct planes. This occurs when the 
Shakura & Sunyaev dimensionless viscosity parameter a is < 0.3 and the initial angle of 
misalignment between the disc and hole is > 45°. The break can be a long-lived feature, 
propagating outwards in the disc on the usual alignment timescale, after which the disc is 
fully co- or counter-aligned with the hole. Such a break in the disc may be significant in 
systems where we know the inclination of the outer accretion disc to the line of sight, such 
as some X-ray binaries: the inner disc, and so any jets, may be noticeably misaligned with 
respect to the orbital plane. 
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1 INTRODUCTION 

In realistic astrophysical situations the angular momentum of a 
spinning black hole may be significantly misaligned with re- 
spect to a surrounding accretion disc. This is expected in black 
hole X-ray binaries where the hole may have received a sig- 
nificant kick during formation jshklovskii 19701: ISutantvo 19781 ; 
lArzoumanian. Chernoff & Cordes 20021 ; iHobbs et al. 20051) . It is 
also likely for accretion on to a supermassive black hole (SMBH) 
in the centre of a galaxy. The scale of the SMBH is so small com- 
pared with the galaxy that any process which drives gas down to 
the galactic centre is unlikely to have any special orientation with 
respect to the SMBH spin. The spin direction is set by the SMBH's 
accretion history and not by whatever causes the next accretion 
event. Processes such as star formation, merging of satellite galax- 
ies, gas cloud collisions and other galactic phenomena likely to 
drive gas inwards are chaotic, so it is reasonable to expect a random 
spread of angles betwe en the angular momentum of infalling gas 
and that of the hole (cf iKing & Pringle 200dKing & Pringle 20071 ; 
iKing. Pringle & Hofmann 20081) . This gas is likely to form a ring 
around the SMBH which then viscously spreads into a misaligned 
accretion disc. 

The Lense-Thirring (hereafter LT) effect of a spinning black 
hole causes non-equatorial orbit s to precess at rates de pendent 
on the distance from the hole l lLense & Thirring. 19181) . In an 
accretion disc, this differential precession causes a warp which 
propagates through the disc. This is the Bardeen-Petterson effect 
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dBardeen & Petterson. 19751) . In this paper we restrict our attention 
to viscous Keplerian accreti on discs with negligible self-gravity. 
IPapaloizou & Pringle (19830 linearised the equations of fluid dy- 
namics to derive the first consistent equation governing the evo- 
lution of a warped disc of this type, assuming a tilt angle smaller 
than the disc angular semi-thickness. They introduced two effec- 
tive kinematic viscosities controlling the angular momentum trans- 
port in the disc. In principle these two quantities depend on the 
nature of the warp: V\ governs the usual radial communication of 
the component of angular momentum perpendicular to the plane 
of the disc due to differential rotation (for a flat disc this is the 
Shakura & Sunyaev alpha viscosity), while v 2 governs the radial 
communication of the component of angular momentum parallel 
to the local orbital plane; t his acts to flatten tilted rings of gas. 
IPapaloizou & Pringle (19831) also showed that conservation laws 
require a relation between vi and v%. For small warps this is 
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where a is the usual dimensionles s viscosity parame- 
ter for discs jshakura & Sunyaev. 19731) assumed <k 1 by 
IPapaloizou & Pringle (19831) . This gives 



a 2 



1 

2a' 



(2) 



(again for small a) where is an equivalent dimensionless viscos- 
ity p arameter for v 2 . 

IPringle (19921) derived an evolution equation valid for larger 
warps expressing global angular momentum conservation for a disc 
composed of rings of gas with no internal degrees of freedom. 
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IPringle (19921) went on to develop a numerical technique for in- 
tegrating this equation over rings of gas to study the viscous evo- 
lution of a time-dependent warped disc. lLodato & Pringle (20060 
used this to study the Bardeen-Petterson effect. They assumed con- 
stant effective viscosities, related as expected for small amplitude 
warps (eq. [TJ. This treatment suggests that a disc warp must al- 
ways straighten itself out on a timescale shorter than the accre- 
tion timescale (assuming only the LT torque and the internal disc 
torques, i.e. no other external torques), because Q implies v 2 » V\ 
for typical values of a. The realignment process enhances the ac - 
cretion rate by inducing extra dissipation dLodato & Pringle. 2o"o6h . 



2 NUMERICAL METHOD 



2.1 Evolution equation 



The analysis of IPringle (19921) did not allow for any inter- 
nal degrees of freedom with in the gas rings, and was d escribed 
as t he 'naive appr oach' in IPapaloizou & Pringle (19831) . How- 
ever lOgilvie (1999b derived equations of motion directly from 
the full three-dimensional fluid dynamical equations with an as- 
sumed isotrop ic viscosity. This confirmed that the equations of 
IPringle (19921) are formally valid only when two differences are 
taken into acc ount (or can be safely neglected). The first is that the 
IPringle (19921) equations omit a torque causing rings to precess if 
they are tilted with respect to their neighbours. In a viscous Kep- 
lerian disc where a <K 1 but not close to zero, this torque is much 
smaller than the usual viscous torques included by IPringle (19921) 
and so can be neglected in time-dependent calculations, (however 
for completeness we include the torque in this work). Our main 
focus, however, is on the second extension. The quantities v\ and 
v 2 are functions of the disc structure and so depend on the warp 
amplitude \if/\ = R \dt/dR\, where t is the unit angular mome ntum 
vector and R is the spherical radius coordinate. lOgilvie (19991) used 
the equations of fluid dynamics to determine the relation between 
the effective viscosities for a general warp amplitude. In contrast 
to an assumption of constant effective viscosities, this analysis sug- 
gests that these quantities drop in a warped region, making it much 
harder for the disc to straighten any twisted rings of gas. Indeed, 
if the warp becomes large enough the disc may break into distinct 
planes with only a tenuous connection. We expect this to occur in 
large warps, as nonlinear cor rections to the fl uid flow are significant 
for warp amplitudes > a dOgilvie. 1999h . 

There have been several theoretical investigations whic h show 
evidence of breaking discs. lLarwood & Papaloizou (19971) report 
an example in Smoothed Particle Hydrodynamics (SPH) simula- 
tions of a circumbinary disc with a strong radial density decrease. 
This is suggestive, since the evolution of a misaligned circumbi- 
nary disc is qualitativel y similar to that of a misaligne d disc around 
a spinning black hole I Nixon. King & Pringle. 201 lh - SPH simu- 
lations bv lLodato & Price (20ld) also found disc breaking. They 
assumed the presence of a large-amplitude warp and followed its 
viscous development. Their results ag reed with those expected for 
the constrained viscosities derived in lOgilvie (19991) . and showed 
that the disc could break if the warp amplitude was large enough 
(note that this paper did not attempt to show that the LT precession 
could produce such a warp). 

In this paper we explore disc evolution under the LT effect 
with general warp amplitud es, and the effective viscosities con- 
strained as in lOgilvie (19991) . We explore the evolution for a range 
of a values and different degrees of misalignment between the disc 
and hole. We aim to discover if the internal torques in an accre- 
tion disc can sustain a smooth warped configuration under the LT 
torque, or whether instead the disc breaks into distinct planes. 



The problem of disc breaking is complex, and to make progress 
we adopt the simplest approach. We use a ID Eulerian ring 
code of the type des cribed in IPringle (1993) and used in 
lLodato & Pringle (20061) . Our code has been adapted to include 
the third effective vi scosity and th e nonlinear fluid dynamical 
constraints found by lOgilvie (19991) . Here 'ID' means that the 
properties a ssigned to each ring dep end only on its radius. Sec- 
tion 2.1 of lLodato & Prin gle (2006) gives a detailed discussion 
of the validity of this approach. The main restriction is that the 
evolution of a disc warp mu st be diffusive rather than wave-like. 
IPapaloizou & Pringle (19831) show that this holds if a > H/R, 
which is usually satisfied for accretion discs around compact ob- 
jects. 

Our ID ring code describes the evolution of warped accretion 
discs by evolving the angular momentum density vector of each 
ring L(R, t) through the equation 
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where vi , v 2 and V3 are the effective viscosities, CI (R) is the local 
azimuthal angular velocity, E (R) is the disc surface density, t (R) 
is the unit angular momentum vector, and il p (R) is the precession 
frequency induced in the disc by the LT effect (defined below). We 
note that L = I,R 2 Slt. 

Equation l[3j shows that there are five independent torques act- 
ing on the rings of gas. The first four terms on the rhs represent the 
internal disc torques responsible for communicating angular mo- 
mentum. The last term represents the LT torque. Here we briefly 
discuss these terms. 

The first term on the rhs of l[3} describes the usual viscous 
diffusion of mass. This is governed by the azimuthal shear viscosity 
vi. The second term is also diffusive. This term, governed by the 
vertical viscosity v 2 , is responsible for diffusing the disc tilt. The 
third term is an advective torque. Depending on its sign this advects 
angular momentum (and hence mass) inwards or outwards through 
the disc. 

The three terms discussed above are exactly those derived 
straightf orwardly from the conservation of mass and angular mo- 
mentum l lPringle. 19921) . The fourth ter m is a precessional torque 
found by the analysis in lOgilvie (19991) . This torque causes rings 
to precess about t x || and thus only acts when the disc is 
tilted. This lead s to a dispersive wave-like propagation of the warp 
dOgilvie. 19991) . The direc tion of the pre cession is dependent on 
a and (cf. Figure 5 in lOgilvie. 19991) . This term has exactly 
the same form as the second term in l[3} with the substitutions 
0.5v2 — » V3 and dl /dR — » ( X dt /dR. This transformation shows 
that we can i ntegrate this term using the technique described in 
IPringle (1992h for integrating the second term. 

The fifth term is the LT torque arising from the gravitomag- 
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Figure 1. The effective viscosity coefficients Q\, Qi and Q3 plotted against warp amplitude \if/\, assuming that a = 0.2, at = and T = 1. The coefficients 
change significantly for < \{fj\ < 5, after which they are approximately constant at a small but non-zero value 



netic interaction of the disc with the spinning hole. This induces 
precessions in the disc orbits, and hence causes the angular mo- 
mentum vector for each ring to precess around the spin vector for 
the black hole with a frequency 
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with c the speed of light and G the gravitational constant. The 
magnitude of J h is given by J h = acM(GM/c 2 ) where a is 
the dimensionless spin p arameter and M is the mass of the hole 
dKumar & Pringle. 19851) . We include the back-reaction on the 
hole angular momentum vector / h by summing the effect of the 
precessions over all radii, giving the torque on the hole as: 
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Since ft p oc J b this equation makes it clear that the magnitude of /i, 
is conserved and hence J b can only precess on a sphere, simply be- 
cause the LT torque has no component in the direction of /h for any 
disc structure. Equivalently the LT torque induces only precessions 
in the disc, so that the reaction back on the hole is simply a sum of 
precessions, which is itself a precession. 

We integrate both <[3j and (0 using a simple forward Euler 
method jPringle. 19921) on a logarithmically spaced grid. However 
we adopt a different implementation of the boundary conditions. 
Rather than implementing a sink over a few grid cells we simply 
remove all of the angular momentum that reaches the inner- and 
outer-most cells at the end of each integration step. This allows an- 
gular momentum (and hence mass) to flow freely across the bound- 
ary, giving the usual S = 0, torq ue-free, accreting boun dary con- 
dition. This is the form adopted in lLodato & Price (20101) . We note 
that this boundary condition does not conserve angular momentum 
exactly. However this small effect does not affect the inclination of 
any of the rings and hence does not alter our results while the warp 
is propagating through the computational domain. Once the warp 
reaches the outer boundary it flows freely off the computational do- 
main and the disc is flattened (i.e. we do not impose a boundary 
condition at the outer edge of the disc to maintain the disc warp). 



2.2 Generating the constrained effective viscosities 

We use the effective viscosities derived bv lQgilvie (19991) for a lo- 
cally isothermal disc. These are constrained to be consistent with 
the internal fluid dynamics of the disc, assuming an isotropic phys- 
ical viscosity. The physics of these constraints is simply local con- 
servation of mass and angular momentum in the disc. 



The effective viscosities which describe the transport of angu- 
lar momentum in a warped disc take the usual form of a Shakura 
& Sunyaev disc viscosity. The constraints determine the nonlinear 
coefficients (gi, Qi, Q3) of the effective viscosities (cf equations 
1 101 1 1 1 l andl 1 2t. These coefficients all decrease in magnitude in the 
presence of a strong warp, where the disc rings are highly inclined 
to each other. Communication of angular momentum is reduced be- 
cause the rings of gas are no longer in perfect contact. We expect 
angular momentum to be communicated through perturbations of 
the particles from circular orbits. The rate of communication, and 
hence the magnitude of the viscosities, therefore depend on how 
easily these particles can interact. This is harder for particles on 
inclined orbits, and once these reach a critical inclination this can 
only happen where the rings cross. We note that for inviscid discs, 
wh ich we do not co nsider here, this interpretation is no longer valid, 
sec lOgilvie (1999T) for further details. 

We now detail the calculation of the effective vi scosities from 
the dis c properties. T he notation in lOgilvie (1999T) differs from 
that in j Pringle (1 992). Our equation [3] and equation 122 from 
lOgilvie (1999f7 relate the effective viscosities (vi, v 2 and vi) to the 
nonlinear viscosity coefficients (Qi, Qi and Q$) by 
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where J* is given by 
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Assuming £ constant in azimuth, and a locally isothermal disc with 
sound speed c s , we can write this as ,JP oc Y,c 2 s /Q. 2 . For a point-mass 
potential with Q. = (GM/R 3 ) [/2 we get after some algebra 
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2Q 2 [(H/R) 2 R 2 n] 



and 

v 3 = Qi [(H/R) 2 R 2 n], 

where H = c s /fi is the disc thickness. 
For simplicity we assume that 
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Figure 2. Disc structures for the constant effective viscosity simulations, given by plotting the tilt angle (0) between the disc and the hole against the log of 
the radius. The initial misalignment of the disc is 8q = 10°. On each plot the legend gives the value of a used in the simulation. From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. Note that in the latter case the warp has moved to the edge of our grid and the 
disc aligns. 




Figure 3. Disc structures for the full effective viscosity simulations, given by plotting the tilt angle between the disc and the hole against the log of the radius. 
The initial misalignment of the disc is 8q = 10°. On each plot the legend gives the value of a used in the simulation. From left to right the panels correspond 
to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. When compared to Fig.[2]this figure shows the negligible effect of the full effective 
viscosities when the misalignment between the disc and hole is small. 



(H/R) 2 R 2 Q. = const. (13) 

This removes the radial dependence of the viscosities, making them 
functions only of \if/\ and a and so allo wing a direct compariso n 
with the constant viscosity simulations of jLodato & Pringle" (20061) . 
The assumption U3\ requires that the sound speed c s oc R~ 3/4 . This 
compares with th e usual steady-state disc, where c s oc R~ 3/i (e.g. 
iFrank et al. 20021) . 

We note that for a di sc without a warp, vy reduc es to exactly 
the usual a disc viscosity dShakura & Sunvaev, 1973T) . since 

v, = ac s H = a{H/R) 2 R 2 £l (14) 

and 

Gi(W = 0) = -3a/2. (15) 

The effective viscosities are now fixed by the nonlinear coeffi- 
cients <2i («> 101), Qi I'/'l) an d Qi (a, |< ft|). We calcula te these co- 
efficients in the same way as described in lOgilvie (1999b withT = 1 
and Ob = using the code kindly provided by Gordon Ogilvie. Fig- 
ure [T] show the coefficients plotted against the warp amplitude for 
a = 0.2. 



from R in = IOR 1SCO to i? D ut = 1000/?i sctJ using 100 grid cells (where 
R lxo = 6GM/c 2 is the radius of the innermost stable circular orbit 
around a Schwarzchild black hole ). The disc surface density is ini- 
tially a Gaussian ring peaked at R = 500R, xa with width 50i? isco . 
We initialise each ring with a constant misalignment angle, do, such 
that the initial disc is flat but tilted to the plane of the black hole. 

lLodato & Pringle (20061) suggest that the only governing pa- 
rameters for simulations of this type are the ratio of disc to hole an- 
gular momenta (Ja/J^), the ratio of the warp radius R w = Cl p R 3 /vi 
to Rq, and the ratio of the effective viscosities V2/V1 (note that they 
neglect v 3 ). However, as we use the constrained effective viscosi- 
ties, V2/V1 is now determined once we have chosen a. 

The parameters that we are free to choose in these simulations 
are a, H/R, Rq, 0q and Ji/Jh- For black h ole accretion discs H/R ~ 
10~ 3 and so we assume this throughout. lLodato & Pringle (20061) 
explored the effect of changing R„/Rq, so we fix Rq as above (note 
that this does not fix the ratio R^/Rq. However this ratio varies little 
throughout our simulations.) To examine disc breaking we choose 
Jd/Jh^ 1 as the simplest case and then vary a and 9q. 



3 SIMULATIONS 

We use the method described above to simulate the warping of a 
misaligned accretion disc under the LT effect. We set up our ini- 
tially planar but misaligned disc on a logarithmically spaced grid 



3.1 Does the disc break? 

We perform several simulations for 0.2 < a < 0.5 to determine 
the various nonlinear responses of the disc to warping. We do not 
simulate a < 0.2 as this is where the viscosity v t is predicted to be- 
come negative in strong warps. We simulate a range of 6q from 10° 
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Figure 4. Disc structures for the constant effective viscosity simulations as in Fig. [5] but this time for initial misalignment 8 = 30°. From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. 




Figure 5. Disc structures for the full effective viscosity simulations as in Fig. [5] but this time for initial misalignment 8 = 30° . From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. When compared to Fig.|4]this figure shows that even for a 30° misalignment 
angle the effect of the full effective viscosities can be small, however there is a steepening of the disc profile at early times for the a = 0.2 simulations. This 
suggests that for the nonlinear effects to be important at 8q = 30° we need a <K 0.2. 



to 60°. This allows us to explore the two regimes: \if/\ < a, where 
nonlinear effects should be negligible, and 1^1 > a, where nonlinear 
effects should be important. Our setup corresponds to J^/J^ ~ 0.02, 
with R w /Rq varying between ~ 0.95 and ~ 1.20 as a changes. Note 
that we evaluate {13} at R with H/R = 10~ 3 . 

We perform these simulations for two cases. In the first we use 
constant effective viscosities. This acts as a control so we can iso- 
late the effect of introducing nonlinear dynamics in the second set 
of simulations, where we use the full effective viscosities. The con- 
stant effective viscosities are simply calculated using the method 
described in Section lZ2l with |^| = 0. 

For the simulations involving constant effective viscosities we 
expect to recover the Bardeen-Petterson effect as seen by sev- 
eral other authors (e.g . iBar deen & Petterson 19751 : IPringle 19921 : 
IScheuer & Feiler 19961 : lLodato & Pringle 20061 etc). However the 
simulations involving the full nonlinear effective viscosities should 
modify this somewhat. 



In Fig.[3]we show the same simulations as in Fig. [2] However 
this time we use the full effective viscosities. We expect the evolu- 
tion here to be similar, as the warp is never able to achieve a large 
amplitude, making the viscosities similar in both cases. This is in- 
deed the case, as there is no visible difference between Figures[2]& 

El 



3.1.2 9 = 30° 

In Figures [4] & [5] we show the disc structures for simulations 
with 6q = 30°. Again there is very little difference between the 
disc structures obtained using the constant effective viscosities and 
those with the full effective viscosities. However the disc profile 
steepens in the early stages of the a = 0.2 simulation. This sug- 
gests that for 9 = 30° we would have to go to a <k 0.2 to get disc 
breaking. 



3.1.1 0=10° 

In Fig.[2]we show the disc structures for simulations with 6 = 10° 
and the constant effective viscosities at 0.01, 0.1 and 1 r v j sc . Here 
we define the viscous timescale as f visc = Rq/vi with v t evaluated 
for \tft\ = and a as defined in the particular case. Fig.[2]shows the 
behaviour for different values of the disc viscosity parameter a. In 
this case the behaviour is similar for all simulations, with lower a 
aligning in a shorter fraction of the viscous time (although in reality 
more slowly because the viscous timescale is longer). 



3.1.3 9 = 45° 

In Figures [6] & [7] we show the disc structures for simulations with 
9 = 45°. There is now a noticeable difference between the constant 
effective viscosity and full effective viscosity simulations. Even for 
large a the disc profile is steepened and for a ~ 0.2 the disc breaks 
into two distinct planes. This break persists for more than a tenth of 
the viscous timescale for the disc. The disc aligns when the break 
propagates out to the edge of our grid which occurs at t = 0.13f v j sc . 
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Figure 6. Disc structures for the constant effective viscosity simulations as in Fig. [5] but this time for initial misalignment 8 = 45°. From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. 




Figure 7. Disc structures for the full effective viscosity simulations as in Fig. [3] but this time for initial misalignment 8 = 45°. From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. When compared to Fig.[6]this figure shows that the nonlinear effects can 
significantly influence the evolution of the disc. Here the profile is noticeably steepened at t = 0.01f v i sc for all a. By t = 0.1? v ; sc the disc break is only 
maintained for small a. For a = 0.2 the break propagates outwards in the disc until reaching the outer edge of the grid (at t = 0,13f v i sc ) at which point the disc 
aligns. For larger a the disc is smoothed before it aligns. 



3.1.4 = 60° 

In Figures [8] & [9] we show the disc structures for simulations with 
9 = 60°. There is again a noticeable difference between the con- 
stant effective viscosity and full effective viscosity simulations. At 
t = 0.01/yisc the disc profile is steepened noticeably more than in the 
9 = 45° case and by t = 0.1? v j sc all of the simulations have main- 
tained a break. We can also see that at t = 0.01/ v isc in the constant 
effective viscosity simulation with a = 0.2 it is possible for the disc 
to break without the nonlinear effects. 

These simulations suggests two main results of using the full 
effective viscosities. First, discs can break for higher values of a, 
and second, discs can break for lower inclination angles. 

There are two features of the simulations which we have not 
discussed so far. The first is the wiggles in the disc at the top of 
the warp in Figures [7] & [9] In these simulations the disc appears 
to try to break into more than just two distinct planes. It may be 
that for smaller a or larger \i//\ the disc can break into more than 
two planes. However we deem this beyond the scope of the present 
investigation. 

Second, the torques on the disc in the strong warp regions gen- 
erate rapid mass flow rates into the ring inside the warp. This results 
in prominent spikes in ~L(R) inside regions with large Increas- 
ing resolution only resolves these spikes better. We therefore leave 
the investigation of this effect to future work, which will require a 
different numerical method. 



4 DISCUSSION 

We have simulated the evolution of a misaligned accretion event on 
to a spinning black hole. Many previous simulations of this type as- 
sumed that the effective viscosities governing the angular momen- 
tum transport were constant. Here we have used effective viscosi- 
ties constrained to be consistent with the internal fluid dynamics of 
the disc. We compare these results with those of equivalent sim- 
ulations which assume constant effective viscosities. In the latter 
case we find that for small enough a the precession induced by the 
LT effect can cause the disc to break into two distinct planes if the 
misalignment is large enoug h. We also find that if we allow non- 
linear effects togilvie. 19991) the disc can break for higher values 
of a and smaller inclination angles. This suggests that nonlinear 
effects significantly change the evolution of strongly warped discs. 
The main effect is in reducing the magnitude of the effective vis- 
cosities in the warped region, noticeably reducing v t more than v 2 . 
This reduces the disc's ability to smooth out the warp induced by 
the LT torque and so prevent a break. This is essentially a modified 
form of the Bardeen-Petterson effect, where the warp is steepened 
by the weakened disc response. 

For small values of a, v% » Vi . Then the effective viscosity 
that tries to flatten tilted rings of gas is much stronger than the effec- 
tive viscosity responsible for transporting angular momentum radi- 
ally. However the vi and V2 torques responsible for communicating 
angular momentum are proportional to t and dt/dR respectively. 
So to break the disc we also need a large warp. As \tf/\ increases the 
V2 torque becomes comparatively stronger that the v\ torque. When 
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Figure 8. Disc structures for the constant effective viscosity simulations as in Fig. [5] but this time for initial misalignment 8 = 60°. From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. 




Figure 9. Disc structures for the full effective viscosity simulations as in Fig. [3] but this time for initial misalignment 8 = 60° . From left to right the panels 
correspond to t = 0.01, 0.1 and 1 viscous times after the start of the calculation. When compared to Fig.[8]fhis figure shows that the nonlinear effects can 
significantly influence the evolution of the disc. Initially the disc profile is significantly steepened for all a and the disc breaks for all a. This break propagates 
outwards in the disc until reaching the outer edge of the grid at which point the disc aligns. 



the discrepancy between the magnitudes of the two torques is large 
enough the disc is effectively unable to communicate radially. As 
the differential precession induced by the LT effect acts faster at 
smaller radii, this allows the vertical effective viscosity (v 2 ) to flat- 
ten each ring in turn into the plane of the spinning black hole. To 
maintain a smooth (small amplitude) warp the disc must be able to 
communicate radially on a timescale comparable to the flattening 
of rings into the plane of the black hole. When this does not hold 
the disc breaks. 

For a ~ 1 the disc can efficiently maintain a smooth warp. 
However for a ~ 0.2 - 0.4 and a large misalignment we show 
there is a strong tendency for the disc to break under an exter- 
nal torque such as the LT torque. We see no reason for this be- 
haviour to change when a < 0.1. However there may be other 
nonlinear effects which become important here. The next step is 
a full 3D hydrodynamical numerical investigation to determine 
the disc response in this regime (e.g. iNelson & Papaloizou 20001 ; 
lLodato & Price 201fj) . 

When the disc breaks it splits into a co- or counter-aligned 
inner disc and a misaligned outer disc. This break can propagate 
outwards and (dependent on a and 6 ) may persist until it reaches 
the outer regions of the disc. At this point the whole disc is co- or 
counter-aligned with the hole spin. This suggests that a break in 
the disc can be a long-lived feature, lasting for roughly the align- 
ment timescale. We note that although we have not included any 
simulations of counter-aligning discs, it is already known that for 
Jd/Jh 1 there is a symmetry in th e b ehaviour of discs of thi s 
type around 6 = 90° jKing et al. 20051 and lLodato & Pringle 20061) . 



Therefore the same behaviour would be seen in discs that are 
counter-aligning with the black hole (i.e. a misalignment of 120° 
would produce the same dynamics with a counter-aligning disc as 
the 60° case produces here). 

For observational purposes it is important to compare the 
alignment timescale f ; ,i with the lifetime fijf c of the system. If 
f ai > hfe» the centre of the disc is never aligned with its warped outer 
plane. In an accreting stellar-mass black-hole binary this would 
mean that jets (which propagate along the axis of the central parts of 
the disc, and so the spin of the black hole) are in general misaligned 
with respect to the or bital plane. This is seen in the black-hole 
binary GRO J1655-40 jMartin et al. 20081 and references therein). 
Simil arly the continuum-fi tting method for estimating black-hole 
spins (Kulkarni e t al. 20lH and references therein) uses the assump- 
tion that the black hole spin is aligned with the orbital axis in esti- 
mating the radiating disc area and hence the ISCO radius. It there- 
fore implicitly requires t a \ < ?ijf c . In our simulations it is clear that 
discs are able to remain warped or broken for significant fractions 
of their viscous timescale. However we note that for strong warps 
there may be o ther nonlinear effects not included in the analysis of 
lOgilvie f 19991) or additional dissipation caused by fluid instabilities 
(cf iGammie. Goodman & Ogilvie 2000l) . These may impact upon 
the alignment timescale. 

It is also possible that the disc is never able to align fully; for 
example an accretion disc in a misaligned black hole binary system 
may be twisted on e way by the LT effe ct and a different way by the 
binary torque (e.g. lMartin et al.. 2009h . We therefore suggest care- 
ful consideration before assuming that black hole accretion discs 
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are aligned - particularly as lLodato & Pringle (20061) suggest that 
the accretion rate through warped discs is significantly enhanced, 
which suggests the possibility that strongly accreting discs are also 
strongly warped. 

Disc breaking can have some important observational conse- 
quences. For example accretion discs in binary systems often have 
an outer disc aligned to the binary plane. However the inner disc 
(and so the direction of any jets) must be aligned with the black 
hole spin, which may be significantly misaligned with respect to 
the binary axis. Similarly breaks in protoplanetary discs may occur 
through the internal disc response to warping during formation of 
the disc, rather than being driven by the presence of a planet. 

In this paper we have used a ID Eulerian ring code to calcu- 
late the evolution of misaligned accretion events on to a spinning 
black hole. We caution that there is still much work to do before we 
can fully understand the evolution of such a disc. For example our 
treatment requires gas to orbit on circular rings, inherently exclud- 
ing any effects which might break the cylindrical symmetry in the 
local disc plane. We also assume that the disc responds viscously, 
with no wave-like behaviour. The method also prevents us from ex- 
amining exactly what happens in the disc break: by its nature, it oc- 
curs over only over a small number of rings and is therefore poorly 
resolved - however this does not affect our conclusion that the disc 
can break. The most frustrating restriction is the difficulty in simu- 
lating small enough a: the simulations we present here suggest that 
for small a nonlinear fluid effects can be particularly important. We 
will address these restrictions in future investigations. 

We remark finally that although we have considered align- 
ment effects induced by the LT effect on an accretion disc 
around a spinning black hole, the torque between a misaligned 
binary and an external accre tion disc has a very similar form 
(N ixon. King & Pringle. 201 lb . Thus we expect similar phenom- 
ena to appear in that case. The enhance d accretion rates through 
the disc due to the dissipation in the warp dLodato & Pringle. 20061) 
would increase the mass flow rate on to the binary. Therefore 
the coalescence tim escale for the binary would be reduced further 
dNixon etal.. 201 lh . 
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